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(57) Abstract: Universal quantum computation is performed using global addressing techniques as applied to a physical system 
of spin entities of two distinct types arranged in an ABABAB sequence. Algorithms are provided for performing global unitary 
operations on such arrays and thereby providing a framework for performing any logical computations. Simple spin | units such 
as free electrons are provided for, as well as molecular entities. The preferred system consists of an ABAB linear array of Group 
V endohedrally doped fullerenes which can be reliably manipulated and arrayed on silicon. Each molecule spin site consists of a 
nuclear spin coupled via a hyperfme interaction to an electron spin. The electron spin of each molecule is in a quartet ground state 
s = 3/2 Neighboring molecular electron spins are coupled via a magnetic dipole interaction. A quantum cellular automata quantum 
computing architecture is provided using these molecules by encoding the quantum information on the nuclear spins while using the 
electron spins as a local bus. The NMR and ESR pulses required to execute the basic cellular automata operations are described. 
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QUANTUM CEIiXiUIiAR AUTOMATA COMPUTING WXTH ENDOHEDRAIi FULIiEREiSnES 

Technical Field 

This invention relates to quantum computing methods, systems and apparatuses. 

5 

Background Art 

Quantum information processing (QIP), relies on the ability to perform a selection of unitary 
operations on a multi-partite system. More precise criteria that must be met by any physical 
implementation of a quantum information processor have been proposed in [1] (see reference 
10 details at the end of the description). In particular, in a spin-based quantum computer 

architecture the spin sites (i) should be easy to physically manipulate, (ii) easy to replicate in 
numbers, (iii) they should possess some type of inter-site interaction and (iv) they also should 
be somewhat isolated from their environment. 

15 Meeting such ideals in practice is very demanding. It is believed by some that a solid-state 
implementation could prove highly advantageous in terms of the scaling up of the number of 
spin sites and many solid-state proposals for a quantum information processor are currently 
under study [3,554]. However, to achieve (iv), usually requires the spin-sites to be located 
within a near-perfect crystal lattice and this leads to great difficulties in satisfying (i) and (ii) 

20 above. 

If the interaction between spin sites is mediated by a predorainantly electronic interaction, as is 
the case in ion-trap quantum computers and [3], the quantum computer operation will be very 
sensitive to stray electric charges. Such systematic charge noise can be source of heating and 
25 may be very difficult to eliminate. 

Additionally the use of local qubit addressing is ubiquitous in very many quantum computer 
designs. However the use of such local qubit addressing comes with a potentially very large 
overhead such as the building of nanoscopic metal gates as in [3], or the greatly increased 
30 bandwidth required in NMR quantum computing in order to fi-equency address a large number 
of spins in a molecule [2]. For small numbers of qubits, this overhead may be acceptable but 
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the decoherence effects of imperfect local addressing in a device containing many tens of 
qubits may be much more serious. 

Finally, the use of nuclear spins alone in a spin-based quantum computer design leads to a 
serious difficulty in achieving a sufficiently high initial spin polarization with which to 
initialise the quantum computer. This is a primary obstacle in liquid state NMR quantum 
computing and the problem remains in solid-state designs even when the nuclei are cooled to 
milli-Kelvin temperatures. 

Faced with all these challenges, the concept of encoding a qubit into the spin system contained 
within an endohedrally doped foUerene has been investigated [6,7]. In [7], the unusual 
properties of Group V C60 endohedral fullerenes are outlined and their potential use in a solid- 
state implementation is indicated. However^ none of the prior art provides a quantum 
computing architecture for use with this material. 

Disclosure of the Invention 

The invention provides a method of quantum computing in a system comprising a plurality of 
spin systems of two distinct types A and B arrayed in an ABABAB. . . sequence, wherein each 
A-type spin system comprises a first spin subsystem Al and a second spin subsystem A2 
coupled to one another, and each B-type spin system comprises a fiLrst spin subsystem Bl and a 
second spin subsystem B2 coupled to one another, wherein the first spin subsystems (AI3I) of 
neighbouring A and B-type spin systems are substantially isolated from one another and the 
second spin subsystems (A2,B2) of neighbouring A and B-type spin systems are coupled to one 
another, said method comprising the steps of 

a) encoding quantum information on the first spin subsystems (AI3I); and 

b) employing the second spin subsystems (A232) as a local bus for the information encoded 
on the first spin subsystems. 

The invention fiirther provides a quantum computer comprising: 

a) a plurality of entities providing quantum spin systems of two distinct types A and B 

arrayed in an AB AB AB . . . sequence, wherein 

each A-type spin system comprises a first spin subsystem Al and a second spin 
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subsystem A2 coupled to one another, and each B-type spin system comprises a first spin 
subsystem Bl and a second spin subsystem B2 coupled to one another, 

the first spin subsystems (Al^Bl) of neighbouring A and B-type spin systems are 
substantially isolated from one another and the second spin subsystems (A232) of 
neighbouring A and B-type spin systems are coupled to one another, 

b) an NMR/ESR pulse generator associated with said array of entities for generating 
and applying said pulses to said array of entities; and 

c) means for reading the spin state of one or more of said entities. 

In another aspect the invention provides a method of quantum computing in a system 
comprising a plurality of spin systems of two distinct types A and B arrayed in an ABABAB 
sequence, wherein each spin system comprises a spin Vz entity, said method comprising 
performing one or more global unitary operations A^ on the A spin systems, wherein: 

A^ is the conditional application of the unitary U when NOT to the A spins 

depending on the state of ^'s neighboring B spins; 

/ is the sum of the states of the neighboring B spins, with spin up(down) defined as 
1(0) such that [0,1,2]; 
and wherein: 

said global unitary operations are performed by applying pulse sequences according 
to one of the following sequence notations (a)-(c); 

a) Ai = H(B, A- {Xf}) . {XV^} . {^f^/J . E(S, A; {X:^}y{Z^^/2} ' {X^/2^s} ' A; {Xf}) , 
read from right to left where 2d = 6; or 

b) = {Z/} - J . {Zf^/J • A; {X^}}{Z^./2} ^ {X^/2^s} ' S(B, A; {Xf}) , 
read from right to left where 5 = 6; or 

c) ^0 implemented as an application of i^^} to all the B spins, followed by followed by 
an application ofW} to all the 2? spms; 
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wherein the controlled operation of applying U to an A spin with left and right B neighbours 

such that the B spins control the unitary application is denoted by C(B^ A; U} and the function 3 
is in turn defined by the equation: 

the left? the right j 

and wherein in equations (a)-(c) the NMR/ESR notation for 7i/2 pulses used is: 

[^S^] ^ exp(i|5f ) , 

and rotation through the angle 7 about the z-axis, for example, is denoted as: 

In this aspect of the invention there is also provided a quantum computer comprising: 

a) a plurality of entities providing quantum spin systems of two distinct types A and B arrayed 
in an ABABAB. . . sequence, wherein each spin system comprises a spin entity copied to its 
two neighbours 

b) an NMR/ESR pulse generator associated with said array of entities for generating and 
applying said pulses to said array of entities; and 

c) means for reading the spin state of one or more of said entities. 

Brief Description of the DrawinRs 

The invention will be further illustrated by the following descriptions of embodiments thereof 
given by way of example only with reference to the accompanying Drawings in which: 

Figure lA depicts the simulated NMR spectra of ^ =^^N@C60 and B =^^P@C60 with (a) 
and (b) labeling the spectral lines from A and B respectively; 

Figure IB depicts the simulated ESR spectra of ^ =^^N@C60 and B =^^P@C60 with (a) 
and (b) labeling the spectral lines from A and B respectively; 

Figure 2 is a schematic depiction of the ABABAB quantum cellular automata chain using 
the Group- V endohedral fuUerenes; 

Figures 3A to 3E collectively form a schematic representation of the steps (a)-(e) 

described below to execute the global operation in a chain similar to that of Fig. 2, having 
the effect of flipping the state of the leftmost A nuclear qubit. 



Figure 3 A shows the S qubits in their ground state with some pattem of / qubits; 
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Figure 3B shows the result of a hyperfine 5Wi4P of the quantum information in the B 
molecules; 

Figure 3C shows the result of the conditional flipping of the A electronic qubit depending 
on the neighboring B electronic qubits; 
5 Figure 3D shows the result of a hyperfine SWAP of the B qubits back into the nuclei; 

Figure 3E shows a controlled global unitary operation U applied to the nuclei of all 
molecules conditioned by the state of their electronic qubits • 

Figure 4A shows an NMR/ESR pulse sequence, in which the solid black bars are 

7r/2 

pulses while the white bars are refocusing n pulses, for a logical SWAP between spin ^ = ^/^ 
10 and the inner electronic qubits of spin = 3/2 described in equation 12 below; 

Figure 4B shows an NMR/ESR pulse sequence, in which the solid black bars are 
pulses while the white bars are refocusing n pulses, for a logical SWAP between spin 
^ = 1/^ and the outer electronic qubits of spin ^ ^ 3/2 described in equation 13 below; and 
Figure 5 shows an NMR/ESR pulse sequence, in which the black pulses are TT pulses 
15 and the, to selectively couple with while decoupling S-^ with S^, and simultaneously 
decoupling S'^ with I A 

Detailed Description of Preferred Embodiments 

20 Quantum Computer Architecture 

Typical quantum computer architectures assume the capability of locally addressing every 
qubit. Such a requirement is highly challenging to engineer and the achievement of this level 
of individual quantum control is a major goal in almost every current implementation. The use 
of such local control methods forces an interaction of the quantum information with very many 

25 "classicaP' control gates, each possibly providing a decohering effect on the quantum 
information. It is also important that the effects of such control gates on the computer's 
Hamiltonian remain as "classical" as possible in that they effect a change in the parameters 
appearing in the system Hamiltonian and do not themselves become entangled with the 
quantum computer. This criteria of classical gating becomes non-trivial as the length scales of 

30 the quantum computer architecture reduces to nanometers. 
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6 

An alternative architecture is to utilize a small number of identifiable spins, placed in a regular 
spatial pattern, and to manipulate the quantum information encoded in this spin chain via global 
addressing techniques. This first of such globally addressed architectures was invented by 
Lloyd in 1993 [26], and utilized three types of addressable spin arranged in the linear pattern 
5 ABCABCABC^ and where each spin site encodes one logical qubit. In fact this cellular 

automata design was one of the very earliest quantum computer architectures proposed in the 
literature. Lloyd showed that such a quantum computer architecture was universal. 

Benjamin [27], found a similar architecture which used only two types of identifiable spin 
10 species arranged in the alternating linear pattern ABABAB. This reduction in spin resources 
came with an increase in logical encoding: a logical qubit is now encoded into four spin sites 
with a buffer space of at least four empty spin sites between each logical qubit. 

The operation of Benjamin's architecture centers on the ability to implement the global unitary 
15 operator, A^. Denoting the spin up(down) state as i(o), Aj is the conditional application of the 

unitary U to the A spins in the alternating spin chain ABABAB, depending on the state of ^'s 
neighboring B spins. In particular, letting / be the sum of the states of the neighboring B spins, 
we have / e [0,1 ^2] Thus, for example, AY is the conditioned application of U to all A spins 

in the alternating chain which have neighboring B spins that are different from each other, i.e. 
20 f = h One has a similar definition for B^, for the application of U on all B spins conditioned on 

their neighboring A spins. The case when U NOT, in Aj (or Bp, occurs quite firequently in 

Benjamin*s architecture and is shortened to Af (or Bf)- 

Benjamin shows in [27], that through a judicious sequence of applications of A^ and B^, one 
25 can implement single qubit operations and two-qubit CNOT operations and thus the 

architecture is universal. In particular, to move qubits through the spin chain one applies an 
alternating pulse sequence of Ai followed by B^, while the generation of a control-U between 
two neighboring logical qubits requires ^ 30 global pulses [28]. One can translate all the 
standard circuit-based quantum algorithms to mn with this quantum cellular automata 
30 architecture [28]. This programming architecture, though somewhat expensive in terms of 
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Spatial and temporal resources, is ideal for use in systems where the individual control of qubits 
is difJicult and it provides a test-bed for various. 

Since Lloyd's and Benjamin's pioneering work little has appeared in the literature with regard 
5 to quantum cellular automata designs for quantum computers [29]. There is described below a 
quantum cellular automata architecture consisting of an alternating linear array of Group V 
endohedrals, ABABAB (schematically depicted in Figure 2), for which the full spin 
Hamiltonian for the pair AB^ (A =^^@C60, B =^^F@C60X can be given as 

10 H/h - gefXeBjf - gifXNBj;^ + A^^S^I^ 

-^JnSfS^ (1) 

15 where I{S), labels nuclear(electronic) spin, and we have made the secular approximation 
dropping terms which do not commute with the electronic Zeeman Hamiltonian since 

The architecture of such a system will be more thoroughly explored below, but essentially it 
20 depends on developing a "quantum algorithm" for generating Benjamin's global operation A^ 

(similarly for gy). We first investigate the possibility of implementing Ji^, using the S 

(electron), spins alone. 

The global operation 

25 In this section we consider an alternating spin 1/2 chain of electronic spins S-^S^S^S^, We 
begin by considering the S-^S^ imit alone where A spin has orJy one neighboring B spin. We 
assume that we know how to perform the controlled operation of applying U to S^, with 
controlling, an operation we denote by C{B, A\ Uy We shall describe below the details 
regarding the sequence of pulses required to execute C(JS, A\ U). We next consider the case 

30 when ^has two neighbors, e.g. S^S^S^. Remembering that one caimot separately address an 
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individual neighboring B spin one can show that instead of perforaiing C{B, A\ U\ the above 
mentioned sequence of pulses now performs 

E(i3, A; U) ~ C(Z3on the left, A] U)C{Bon the right, A] U) , (2) 

5 where Eq. (2) defines the function 3. One can observe that the case when the desired global 
operation is Ai = A^^^, we have Ai = E(5, A; NOT). 

We now indicate the procedure for implementing when the desired U ^ exp(i03^\ For ease 
of notation we define the standard NMR/ESR symbol for 7r/2 pulses [30531532], 

and the more general rotation through the angle 7 about the axis, for example, by the symbol 

{Z:^} = eiqp{ijSf) . (4) 

We find that we can execute A^, when U ~ ex.p{i9S^) through the following pulse sequence 
1 5 (read right to left), 

= 3(8, A; {Xf}) ■ {X;^/^.s} ■ {^ll/a} ■ 2(^, >t; {A';:^}) 

■{Z^^/2}-{X^/2-s}-SiB,A;{Xf}) . (5) 

20 where 25 = 6, To implement A^ one sets S ~ 0 and replaces the left most term in (5), 

B(5, A; {Xf}), by {Z;^}, while to implement A!q, one applies {X^} to all the B spins, then A^, 
and finally flip the B spins back with {X^^y The above pulse sequence (5), can be easily 
modified to suit other forms of U besides exp{i0S^), however we will find below that the 
most useful case is when ^ = tt and U = exp{i7rS^) NOT on the >t spin. 



25 



All EiectroBic Quantum Cellular Automata Quantum Computer 
From Figure 1, we can separately address the electronic spins S-^ and with a 
multifrequency selective soft pulse of length > 20ns, If we first consider both S-^ and to 
both be spin-1 /2y then one can use standard pulse encodings of the CNOT [3132], 



30 
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9 

C{A, NOT) = ['-S^][:^Sf T S?^][±2SfS^,][-^S^]e^'-/^ , 

C[B, A NOT) = [S^M^Sf T S^][±2SfS^][+S^]e^-l^ , 
5 to generate the operation Ai and i5i of the previous section. 

To perform universal quantum computation one must be at least capable of implementing a 
two-qubit CNOT and a minimal set of one-qubit operations such as the phase gate, the tt/S gate, 
and the Hadamard gate. In the quantum cellular automata architecture of Benjamin this means 
10 that one must at least be able to implement the global operation where U is taken from 

this minimal set of one qubit gates and the case where U ^ NOT. NMR and ESR pulses 
effectively generate terms in the Hamiltonian which are linear in the spin operators, Ia:^y^z^ and 
Soc,y,z' particular, the spin operators when commuted close to form the Pauli algebra. In the 
case of spin- 1/2, the Pauli algebra generates the group SU(2), or all possible spin- 1/2 unitaiies. 
15 Thus one can generate any desired spin- 1/2 unitary through a sequence of pulses when acting 
on an isolated spin. Furthermore, in the case of the alternating chain ABABAB, one can find 
suitable alterations of (7), when used in (2), and (5% to generate any given unitary U in Aj 
and similarly for Bf* Thus an all-electronic quantum cellular automata quantum computer 
architecture (QCAQCA), using a spin-i /2 ABABAB chain is possible. 

20 

Pulsed ESR machines of the type required to generate the required pulse sequences discussed 
above are of course well known in the art. The spin i4 chains can be implemented as defect 
centres in crystals (e.g. diamond) or as single electrons trapped in electrostatic traps in ultra- 
high vacuiim. Further embodiments will suggest themselves to the skilled person and new 
25 implementations will undoubtedly become available as research in this field progresses. 

Endohedral Fullereiie-based Quantum Computing 

In an altemative embodiment, the invention provides a quantum cellular automata quantum 
computing architecture (QCAQCA), which can successfully be applied to operate with Group 
30 V C60 endohedrals. The QCAQCA design only requires an homogeneous constant magnetic 
field. More importantly, in the QCAQCA^ fiiU scale universal quantum computation with 
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10 

quantum error correction only requires frequency addressing at a fixed, small number of 
frequencies, independent of the size of the device. These two advantages mean that a quantum 
cellular automata based device may be far easier to develop in hardware than any locally 
addressed architecture. 

5 

For the design developed below we exliibit the specific pulse primitives necessary for quantum 
cellular automata quantum computation. We are also able to estimate rough figures of merits or 
number of quantum gates that can be executed in a decoherence time. The estimated figure of 
merit obtained below compares well with other proposed quantum computer designs. Again, 
10 pulsed NMR/ESR machines of the t5^e required to apply the necessary pulse sequences to the 
array of endohedrally doped fuUerenes are well known. An example of such a machine is the 
Bmker Elexsys ESR machine. 

Below we review the properties of »-Y@C60 (Af — ^-^N, ^^P), pertinent for quantum infomiation 
15 processing. We briefly outline the essential ingredients of QCAQCA. We develop a number of 
tools (pulse sequences), for manipulating the Hamiltonian of an alternating linear chain, 
ABABAB, of th^e endohedrals where A =^^@C60 and B ==?^P@C60. We consider an 
architecture which also includes the nuclear spins and show, by developing new gate pulse 
sequences, that universal quantum computation can be achieved via a quantum cellular 
20 automata design. We finally, propose a number of readout techniques that are applicable to the . 
proposed ensemble (a T^^e-II quantum computer design [9]), and single-issue (or a Type-I 
quantum computer design [9])^ quantum computer design. 

Group V Endohedrals 

25 The endohedrals ^^@C60, ^^P@C60 and N@C70 all exhibit very sharp ESR spectra [10,1 1]. 
This indicates the presence of free electrons within these molecules. It has been further shown, 
both theoretically and experimentally [12,13,14], that the trapped atom ( iV or P% sits at the 
geometric center of the fuUerene cage and that the electrons are in a 5 — 3/2 quartet ground 
state. 



The trapped atoms are extreme examples of compressed atoms [15^12], where the electronic 
wavefunction of the trapped atom is repulsed away from the encompassing carbon cage and 
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suffers a spatial compression. This unique type of '"bonding"' leads to significantly higher 
nuclear-electron waveflanction overlap and thus larger Hyperfine coupling constants than what 
is found in ''firee" atoms. For the case of C60 Group V endohedrals, the large symmetry 
renders this hyperfine coupling highly isotropic. Another consequence of this compression is 
the virtual lack of any electronic interaction between the trapped atom and the carbon cage. 
The trapped atom is motionally confined by a harmonic-like potential to the center of the cage. 
In almost all respects, the trapped atom behaves as a "j&ee" (unbonded) atom, though spatially 
restricted to be within the fuUerene cage. All of the above findings imply that such endohedrals 
behave as nanoscopic molecular neutral atom traps. 

The distributed tt bonding electrons on the C60 also act as an almost perfect Faraday cage^ 
strongly isolating the electrons of the trapped atom firom external electric fields [16]. This latter 
observation implies that it would require the application intense local electric field gradients 

IV /nm, in order to alter the Hyperfine coupling constant of the trapped atom. Such an 
electrostatic addressing scheme has been proposed to address nuclear spin qubits in the 
Phosphorus donors of the Kane design [3], but with far lower field gradients. Such a scheme 
would prove very difficult to execute here due to the requirement of such intense electric field 
gradients. 

One very significant advantage of containing the spin site inside a relatively large C60 
molecule (diameter ^ Inm), is that such endohedrals can be nano-positioned using current 
scanning tunneling microscope techniques, and moreover neighboring C60 molecules can be 
nano-positioned on the Silicon(100)-2x surface to be as close as ^ l.lnm. The C60 molecules, 
on the Silicon(100)-2x siirface, are fixed translationally at room temperature and also 
rotationally at lower temperatures. Due to the electronic- wavefunction spatial compression 
experienced by these molecules, the chemistry of the doped Group V fullerene material is 
practically identical to undoped fullerene. It is therefore possible to self-assemble large 
organised spin structures using existing well-studied synthetic supra-chemistry techniques [7]. 
An altemative implementation is to use endohedrally doped fuUerenes positioned in nanotubes. 

The very sharp ESR spectra fi*om these molecules indicates very long relaxation times and 
more recent measurements have shown the electronic relaxation times are Ti ~ Is at 
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T 7° if, while T2 ^ 20jus, VT [1 8]. No nuclear relaxation times have yet been recorded 
but they are expected to be several orders of magnitude longer than the electronic relaxation 
times. The measured T2 time contains contributions from several sources, e.g. unresolved 
dipolar couplings. The theoretical maximum of T2, in the complete absence of the unwanted 
interactions with other paramagnetic impurities, is the relaxation time Ti- Indeed, for the case of 
Phosphorus defects in isotopically ultra-pure ^Si, the phase relaxation time of the loosely 
boxmd Phosphorus electrons can be as long as r2 > .1 msec with Ti > 1 hour [19], There T2 is 
limited by Hj^erfine interactions with residual ^^Si nuclei [20]. The indicated increase of T2 
towards Ti in the situation of vanishing spin density [22], is crucial for all current spin-based 
quantum computer proposals [21,35455,6,7,23]. Achieving such limits will be very challenging. 
Making these assumptions for the endohedral electrons we expect T2 isec, at 

T T'^K, with perhaps longer times at lower temperatures. 

Such electronic Ti relaxation times (r^^ isec), in such complicated molecules is highly imusual 
and may be unique in all non-crystal hosted spin-sites. Long relaxation times also add to the 
usefulness of such material as hosts for storing and manipulating quantum information. 

An essential ingredient for any quantum information processing is an inter-qubit interaction 
which can generate entanglement. For Group V C60 endohedrals this cannot be a direct 
electronic exchange type interaction as the electronic wavefiinctions are tightly compressed to 
be within the C60 cage. The C60's Faraday cage does not restrict magnetic interactions and 
neighboring endohedrals experience a significant magnetic dipole coupling, 
Hn /(I ~ 3 cos^ d)[SI^S^ - I • S] I^S^, in the weak coupling limit. The strength of this 
dipolar coupling has been estimated to be J^? ^ 50MHz(l/r^), where r is the separation between 
the neighboring trapped atoms in the endohedrals measured in nanometers [24], and 0 is the 
angle between f and the external B field. The Hyperfine coupling constants for ^^@C60, 
^^P@C60 have been measured [25], and are given in Table 1. 
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Table 1 : Table of Hyperfine coupling constants and interaction energies for interacting ^^@C60, 

3ip@C60 molecules at = 2T. 



From the above, to a first approximation, the full spin Hamiltonian for the pair AB, 
(A =^^@C605 B =^^P@C60)5 can be given by equation (1) which is repeated here: 



10 



H/h = ge^eBjf - giiiNBj^ + A^Sfif 



4- JuSfSf (1) 

where I{S\ labels nuclear(electronic) spin, and we have made the secular approximation 
15 dropping terms which do not commute with the electronic Zeeman Hamiltonian since 

^Lrmour - 9efieB, » ( Jd, ojia^rmour - QnP'nB,, A). The individual ESR and NMR spectra of 
these molecules. are simulated in Figure 1, and have been measured in [1 1]. 



Quantum Cellular Automata 

20 The endohedrals ^^@C60 and ^^P@C60, unlike the all-electronic spin Vi chains for which the 
global unitary operators were developed above, possess electron spin 5 = 3/2 and the 
application of Benjamin's QCAQCA is not as straightforward. 

One possibility is to consider the four electronic levels of the electron spin with 5 = 3/2 to 
25 constitute a qu-dit with rf = 4. Multi-level quantum computation has been proposed elsewhere 
[33], but this invention does not extend Benjamin's formulation of the spin-1/2 QCAQCA to 
operate with qu-dits. Indeed, as shown below there are reasons to believe that that such an 

extension may not be possible within the confines of NMR and ESR. 

30 Remaining with qubits we instead propose to encode a qubit into the four levels of a spin 3/2 in 
two ways: an "iimer qubit" {\ms = ±1/2)}, and an "outer qubit" {\ms = ±3/2)}. By a small 
alteration of (7), one can obtain 
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Cir.nAS^. S^; NOT) - C^teriSA. ^ NOT) - [S^][-Sf ^ Sf][+2SfS^][--S^]e+'^/^ , 

(8) 

5 

CinneriS^. S^l NOT) - C^,er (5^, S^; NOT) - [S^]hSf - 5f ][+25/5f ][-5^]e^-/^ , 

(9) 

where dnner is a CNOT on the Hilbert space spanned by the inner qubits, and similarly for 

Couter^ 5^ = NOT{S-'% and we have used the spin-3/2 representations of the SU{2) group. 
Armed with tliis one can easily construct the Benjamin global operations Ai and Bi on the inner 
qubit subspace. One can further construct (as discussed below), the global operations A^^, 
(and similarly Bf^'^)^ for / = 0, 1 and / = 2, the operations most frequently used in the quantum 
cellular automata quantum computation. 

15 

However one finds that one cannot generate all the one-qubit gates required for universal 
quantum computation. As above, ESR pulses essentially generate terms in the Hamiltonian 
that are linear in the spin operators and again these terms commute to form a Pauli algebra. 
Evolution under such terms lead to dynamics that resides in the group generated by this Pauli 

20 algebra (or SU (2)). However, since now we have S = 3/2, this SU{2) is merely a subgroup of the 
full group of unitaries that can act on a spin-3/2^ but as long as the system Hamiltonian is linear 
in the spin operators we need only to consider the dynamics within the SU{2) subgroup. 
Noraially a unitary operation U e SU{2) c SU{4) will not factor into unitary operations on the 
inner and outer qubit subspaces. However if U is diagonal in the \ms) basis it can so factor 

25 and with this one can find pulse sequences to implement where U is a phase gate, and the 
tt/S gate but cmcially not the Hadamard gate. Thus an all electronic-QCAQCA with a system 
Hamiltonian which is linear in the spin operators and where S > 1/2 seems not possible. 

One way of implementing SU{2) operations on the qubit subspaces and one which is well 
30 known in NMR and ESR is to introduce terms in the system Hamiltonian which are not linear 
in the spin operators. Such terms lie in the full spin-3/2 algebra and the generated unitaries no 



wo 03/081526 PCT/EP03/02997 

15 



longer lie in the SU{2) snbgi'oup but in the much larger SU{4) group. One effect of such 
nonlinear tenns, for example Zero-Field-Splitting tenns: Hzfs oc can be to lift the 
degeneracy of the microwave transition frequencies between the \ms) ^ \ms + 1). This allows 
one to frequency address the transition jm^ = -1/2) ^ \ms = +1/2), or the inner qubit subspace, 
5 and directly implement SU(2) operations on this subspace. This allows the generation of 

arbitrary one-qubit unitaries t/ [34]. However, this solution frustrates the coherence transfer 
portion, [+2SfS?l of the pulse sequence (8), and (9), and thus the construction of the global 
QCAQCA operator A^. This can be seen by realising that to obtain a [+2S/5f ] pulse one must 
effectively cancel out all the terms in the system Hamiltonian bar the term SfS^- As standard in 

10 NMR and ESR one does this cancellation through "average Hamiltonian" theory [35], via the 
application of various SU{2) unitaries by RF or MW pulses. As long as the system evolution 
remains within the SU{2) subspace such an averaging out can be achieved. However when 
terms nonlinear in the spin operators are introduced into the system Hamiltonian no application 
of 5(7(2) unitaries can, in general, average out die resulting "nonlinear" evolution (i.e. in the 

15 average Hamiltonian theory nonlinear terms in the Hamiltonian will almost always average to 
other nonlinear terms), and thus the generation of{^2SfSf\ and from this A^, is not possible. 

Thus we have found that when operating on an ABABAB, spin-3/2 chain with ESR pulses one 
cannot generate all of the necessary global operations required for universal quantum 

20 computation. The present invention, however, demonstrates that the combination of nuclear and 
electronic spins witliin the chain can provide the necessary framework for the global 
operations. 

Nuclear-electron quantum cellular automata 

25 ^^@C60 and ^^P@C60 both have nuclear spin-l /2. In this section we assume that the quantum 
information is stored in the nuclear spins of the dopant atoms in the ABABA chain. This has 
the significant advantage that the nuclear spin relaxation times are typically longer than the 
electronic spins by several orders of magnitude. As the nuclei only are coupled via the 

Hyperfine interaction to the electrons, one can, via a suitable RF and MW pulse sequence, 
30 cancel out this Hyperfine interaction with high precision. This allows the nuclei to act as a 
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quantum memory and store the quantum infomiation for very long periods of time in between 
processing. 

15 

To execute the Benjamin global operation Jjj we essentially use the inner electronic qubits as a 

5 local "bus" for the quantum information stored in the nuclei. The procedure can be thought of 
as a quantum algorithm for the generation of the global operation and consists of several steps 
which are pictorially represented in Figure 3. 

The algorithm starts with the inner electronic qubit bus in the ground state (|m^) = | — 1/2) for 
10 all molecules), and this bus is returned to the ground state after implementing the global 

operation JJj. We first (a) assume that we have an arbitrary pattern of quantum information 

written onto the nuclear spins of the A and B molecules with the electrons set as above. We 
then, (b) swap the quantum information of all the B molecules from their nuclei to their inner 
electronic qubits. We then (c) tag those A molecules which will receive the unitary operation 
15 ?7 in Ji^ by flipping the state of the electronic qubit of all the A molecules conditioned on the 

state of it's neighboring B electronic qubits. We then (d) undo (b) by swapping the quantum 
information back into the B nuclei from their electrons. 

At this point all the inner electronic qubits are back in their ground states except for the 
20 "tagged" A molecules. We then^ (e) perform a controlled- U operation on the nuclear qubits 
of all the molecules using their inner electronic qubits as the control. Finally (f), we undo the 
operations (d)^ (c) and (b). The system is then ready for the execution of the next global 
operation. 

25 One restriction we make in the following is that we rule out the use of selective pulses on the 
nuclear spins of the Nitrogen and Phosphorus. From Figure 1, to execute such selective pulses 
would require RF pulses of veiy long duration to frequency differentiate between the NMR 
transitions of the Nitrogen and Phosphorus, Such slow nuclear pulse sequences would lead to 
very long gates times. We will therefore not allow ourselves to use such selective nuclear pulse 

30 sequences and instead make use of fast "hard" nuclear pulse sequences [36]. We do allow 

selective addressing of the electronic transitions of the Nitrogen and Phosphoras as these pulse 
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sequences are much shorter in duration. With this restriction however we are still able to 
design all the pulse sequences that are required to carry out the steps (b)-(e), above. The 
resulting combined total pulse sequence is somewhat lengthy but could be substantially 
shortened by applying the principles of optimal control theory [37]. 

We now expand on the above steps (a)-(e), in detail. In section (A) below we describe the 
initialization of the spin-chain. In section (B) we design &e SWAP pulse sequence required for 
step (b). In section (C) we design the electronic pulse sequence to implement the conditional 
operation required in step (c), and finally in section (D) we outline the execution of the 
controlled- U needed in step (d). In section (E) we discuss the undoing of the previous steps 
(b)-(d), while in section (F) we provide some estimates of the pulse durations and logical QCA 
gate times. 

A:- Initializatioii 

The QCA architecture may be implemented on an ensemble (Type-II) or single-issue (Type-I), 
quantum computer. The achievement of high nuclear polarizations has proved to be a major 
obstacle in NMR quantum computation. It remains a source of difficulty here as well. 
Polarizing the electrons is far easier, and one can achieve a difference in the ground to excited 
electronic populations of 6 0.999, at a temperature of rK and B lOT. Even when multiple 
spins are tensored to produce a pseudo-pure state, with such high individual polarizations one 
can still achieve a pseudo-pure state with a purity on the order of e .998, with 1000 spins at 
these temperatures and magnetic fields. This polarization can be transferred to the nuclei via an 
INEPT pulse sequence [38]. To repump the electronic polarization and achieve a subset of 
completely spin (nuclear and electron), polarized molecules one could consider two chains of 
equal length: an ABAB • ■ ■ chain and a CVCV - • * chain where all A, B, C, and P are globally 
addressable. Since at l^K and S = lOT, half of the spins (the electrons), are completely 
polarised, one can consider an ABAB - - - ABABXCVCT) • - • CVCT^ super-chain where X has a 
different resonant frequency yet again. One can then use a spin cooling quantum algorithm 
[39]^ to efficiently spin cool all the nuclear and electronic spins on one side of X (i.e. X 
acts as a spin gate). One would then switch off the interaction between the two sub-chains at 

X and use the spin polarised half chain to perform the quantum computation. Such a design 
would be suitable for a Type-II (ensemble), quantum computer. Altematively, in the case of a 
Type-I quantum computer one can replace the CVCT> chain with an ABAB chain and place the 
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readouts on the polarised half chain. More practically any Type-I quantum computer will 
require a single-spin readout of the electronic spins of the endohedral atoms. Such a readout 
may also serve to initialize the electronic qubits. We discuss possibilities for a single-spin 
readout further below. 

5 

B:- Nuclear-Electronic SimP 

A cmcial ingredient in the above mentioned quantum algorithm to generate is the Hyp^ne 

SWAP. This operation performs a logical swap between a qubit stored in the nuclear spin and a 
qubit stored in the electron spin. Below we find two types of Hyperfine SWAP, one that 
10 performs the swap using the ixmer electronic qubit and one that uses the outer electronic qubit. 
For the most part we will use the former type of swap but we will find later that the latter swap, 
which incorporates the outer electronic qubit, may be of significant use in the problem of qubit 
readout and we discuss this below. 

15 To develop the pulse sequence for tlie Hyperfine SWAP it is useful to break down the SWAP 
operation into three CNOT operations. Denoting SWAP {I : S), we have^ 

SWAP{I : S) - CNOT{I', S)CNOT{S; r)CNOT{I; S), where CNOT(I; S) is the controlled NOT 
for I controlling 5. In our case 1^1/2 while S = 3/2. Substituting the spin-3 /2 
representations for S into the standard spin-l /2 CNOT pulse sequence (7), fails to yield a 
20 proper CNOT operation. One instead must make reference to the two cases of a CNOT with 
respect to the inner and outer qubits of S, With this one can find: 



25 



30 



CNOT^ter{I\ S) = lSy][:^S, =F Q[±2SJ,][^Sy]e^'-/^ , 



CNOT^uriS; I) - [Iy][±S, ± J.][±2S,/,][-^/^]e^-/^ , (10) 



while 



CNOTtnnerill S) - [Sy][::fS, T Iz][±2SJ,][-Sy]e^'^^^- , 
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With these one can derive 

(12) 

SWAPouteril : S) - [-Iy][-h2I,S,][^Iy - Sy][~I, ~ S ,][-^2I ,S ,][S ^ + /.][+2I,5.]h/^]e+^^"/^ , 

(13) 

and the related pulse sequences are graphically shown in Figure 4. We note for later the curious 
case of (13), where we have in some sense "amplified" the magnetic signature of the quhit from 
a Amj = ±1, transition to a Ams = ±3, transition. In the following we will restrict ourselves to 
the inner electronic qubit subspace and thus (12). 

Having found a SWAP operation between the nuclear and inner electronic qubits of either 
molecule we now must show how one can adapt the sequence (12), to swap out only the B 
qubits while leaving the A nuclear qubits alone. We now consider the case of two coupled AB 
molecules with the Hamiltonian (1). The desired selective swap action can be achieved hy 
replacing the [+2hS.] terms in (12), with the pulse sequence shown in Figure 5. This new 
pulse sequence effectively averages out the Hyperfine coupling interaction in the A molecules 
thus tuming off the action of the Hyperfine CNOT and SWAP gates for these molecules. The 
Hyperfine interaction for the B molecules is left intact by the pulse sequence but is reduced in 
magnitude thus lengthening their Hyperfine gate durations. 

The black pulses in Fig. 5 are TT pulses and the duration of the entire sequence is matched to 

in (13). The gray pulses are fast hard ^tt pulses which are required to extend the 
duration of the [2lzS^ for the case when s wrapping between the nuclear and electronic spin in 
^^F@C60. By placing these pulses at the points - l)/(4m)^ and + l)/(4m)^ in each third 
section as shown, we can extend the pulse sequence by a factor m. Setting m = 6 brings the 
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^^P@C60 SWAP to be ^ 2OMH25 and thus long enough for the electronic spin selective pulses 
to operate on S-^ and S^. 

This particular pulse sequence has a number of advantages, the most important of which is that 
5 we do not selectively address the nuclear spins and thus can apply fast hard pulses to and 

simultaneously. We only take advantage of the selective addressing of the electronic spins of A 
and B as outlined previously. Altematively, one can use soft multi-frequency nuclear pulses in 
the following to approximate the hard nuclear pulse. In more detail, from [40,41], the action of 
this pulse sequence is to remove to all orders (within the secular approximation made 

10 previously), all interactions between and 5*^ ^ J*^. This sequence also removes all the 

Zeeman terms in the AB Hamiltonian and one is thus left with H/h A^S^If. When inserted 
into (12), the effect is to swap the qubit that was stored in 1^, into the inner electronic qubit of 

S^. As mentioned above, the quantum information stored in is not swapped into S-^, but 
these spins do receive the local unitary transformation, The local unitary [5^] is 

15 removed with a spin selective electronic pulse [S^l however we cannot remove the [/^] 

without performing a spin selective nuclear pulse, an operation we wish to avoid. We will see 
however, that this extraneous local operation on /-^ is not important for the tagging of the A 
molecules and furthermore the operation will be "undone" in step (d) (c.f. Figure 3). 

20 C;- Electronic tagging of / 

Summarizing the previous steps: we have now swapped out the nuclear qubits of the B 
molecules in the ABABABAB chain into their inuner electronic qubits. Additionally, all the 
inner electronic qubits of the A molecules are in their ground states jm^) = | - 1/2), i.e. we are 
at step (b) in Figure 3. We will now electronically "tag" those A molecules targeted by the 
25 global operation Af* 

More precisely, we wish to flip the state of the A inner electronic qubit subject to the function 
/, with the control being A's nearest neighbor B's inner electronic qubits. We noted at the end 
of the section entitled "quantum cellular automata", that many inter-molecular electronic 
30 unitaries cannot be achieved with 5*1/(2) pulses. However, the case of / = 1 is simply the 
operation H(B, A; NOT), where one uses (8), for the CNOT in the inner electronic qubit 
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subspace. The case / = 2 is more difficult and one must make foil use of the construction given 
in (5), to obtain the following pulse sequence 

S(B, ^; ^/NOf) - J • 3{B, A^ VnOT) 

where H(i3, .4; \/]VOT) = 0(3^, S^; E)^, (see (2)), and 

CiS^, S^; S) = e4-^>^e-^^^^^+^i-^-"^^§>e-^l-^^^^e^'-/^ . (15) 

10 

With this pulse sequence one achieves the required spin flip of the A molecule's inner 
electronic qubit up to a phase. This global phase is unimportant as it will be removed in step 
(e), when we undo the unitaries. Finally the case / = 0 is dealt with by first flipping all B inner 
electronic qubits, then using the f = 2 operation above, and then flipping back. In the case of 
15 Figure 3, we have chosen to implement the tag / = 1 operation in step (c). Once we have set the 
electronic tag on the appropriate A molecules we "undo" the SWAP, returning B's quantum 
information back into its nuclei and B's iimer electronic qubits to the gi^ound state, arriving at 
step (d) in Figure 3. 

20 D:- Implementing U 

We now can implement any desired unitary on the flagged A qubits. The construction of this 
operation is well known [42], and consists of applying hard pulses [36], to all of the A and B 
and nuclei interspersed by CNOT gates between the inner electronic qubits and nuclear qubits 

of all the molecules. More precisely, the controlled execution of any given unitary operation 
25 on the nuclear spin / controlled by the inner electronic qubit S on all the molecules may be 
written as 

I;U)=^D Oil, S; X) E C{I, S;X)F , (16) 



30 



U^e'''R,{j3)Ry{j)RAS) 
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^e^'^DXEXF , (17) 

In (16). /3, 7, and 6 are the Euler angles of the desired nuclear SU{2) operation, C(Iy 5; X) is the 
Hyper fine control-NOT given in (1 1), the nuclear rotation operators are given by 
5 Rsid) = ejcp{-i6 5/2\ where S =^x,y,z and standard pulse sequences for these spin rotations can 
be found in [38]. Further D = R,(/3) Ry ('y/2), E = Ry{'-j/2)R^(-(S + /3)/2), F ^ R,((S ~ /?)/2), 
X ^ NOT, and a is a phase. Crucially one can show that DEF = I. From this one can see 
that the action of the sequence (16)5 on the nuclear spin is unity when the inner electronic qubit 
is not set. When the inner electronic qubit is set however, the action of (16), is to execute U 

10 from (1 7), on the nuclear spin up to the phase a. One then finally applies the local unitary 

eyip(j.a{Sf -f S^)\ to all molecules to yield the final phase factor a in (17). With this one has the 
very powerful capabilit>^ of applying any desired single-qubit unitary to the nuclear qubits of 
the tagged A molecules. However^ the typical operation in the QCAQCA is Af, and this can be 
executed more simply by a Hyperfine CNOT between the electronic and nuclear qubits of all 

15 the molecules. In Figure 3, we illustrate the application of Ax in subfigure (e). 

E:- Undoing the unitaries 

Following the arguments made above, we execute step (e), that is: once we have implemented 
the desired unitary on the A nuclear qubits we ''undo" steps (b), (c) and (d). We SWAP the 
20 nuclear qubit into the electrons on all B molecules, undo the tagging of the A molecules, and 
then SWAP the quantum information back into the nuclei of the B molecules. The system is 
left with all the inner electronic qubits in their ground states while the nuclear qubits have 
received the global operation JPj. 

25 F:- Gate duration 

As we noted above, the resulting pulse sequence is quite lengthy. However it may be possible 
to compress man}^ of the above operations. Also all nuclear pulses are fast hard pulses of little 
duration [43], while the primary slowdown arises Jfrom the nuclear-electronic SWAP whose 
duration is limited by the value of the Hyperfine coupling constant of ^^@C60 ( ^ 20MHz), 
30 and the separation of the spectral lines of the A and B molecules in the ESR spectrum of Fig. 1 
(-- 50MHz). The typical QCAQCA global operation, Af, entails a pulse sequence in which the 
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terms [+a^b:,\ occurs 15 times, each with an average duration ofbOns, thus roughly bringing 
the cycle time of the global operation down to l/j,s. As noted above, the simplest quantum logic 
gate in the QCAQCA scheme requires ^ 30 global operations. This finally brings the cycle 
time for logical gates in the resulting QCAQCA to be roughly 30/i5. 

As noted above, one has Ti Is at r = 7°iir, and 20/iS for concentrated samples of 
^^@C60 Group-V endohedral material. As mentioned above, recent experiments have 
indicated a dependence on concentration for T2 and it is commonly expected in all cuiTent spin- 
based quantum computer implementations that T2 will rise towards the value of Ti in the limit of 
zero spin concentration. If r2 ^ is then one can expect on average 10, 000 - 30^ 000, logical 
operations within this dephasing time (figure of merit). Indeed, at temperatures lower than T'^K 
one might be able to achieve far greater figures of merit as the decoherence and dephasing 
times increase. A figure of merit of 10^ compares very favorably with most alternative 
proposed implementations for quantum information processors. 

Readout 

One of the most challenging aspects of any solid-state implementation is that of readout. We 
first discuss the possibilities for an ensemble readout for the QCAQCA and later a single-spin 
readout for a Type-I implementation. 

In the original scheme of Lloyd [26], the only spin sites which could be individually addressed 
were those at either end of the chain as these spins were firequency differentiated by having a 
single neighbor. This is also the situation here. However, if these two sites were the only 
readout sites on the chain then the architecture would not be scalable in the presence of noise. 
One ensemble readout possibility is to attach an electrically isolated paramagnetic adduct to 
every B molecule, ( Gd@C82 may be a possibility). Such readout sites would typically possess 
short relaxation times and one would have to isolate such sites from the operation of the 
primary processor using ESR pulses. It should be possible within the cellular automata 
architecture to selectively readout the state of a single qubit, that is, to transfer the state of a 
given logical qubit into the readout site. This can be done by adapting the single-qubit unitary 
operator U, pulse scheme in [27]. This pulse sequence is arranged so that a logical qubit 
(which is usually encoded in four spins, ABAB), gets driven into the state of a single spin 
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which is then subject to a unitary operation, through the application of the global 

operation (say B^, when the qubit has been forced into the spin state of a S molecule). The spin 

neighborhood of this target B spin is arranged so that this global operation only effects the state 
of that spin alone with an appropriate choice for /. For readout one could instead apply an 
operation Bj^'^^, which swaps out the state of the spin into the readout site. This would be read 

and then SM-^^Ped back. 

For a Type-I quantum information processor based on the ABAB chain of ^^@C60 and 
^^P@C60 one must be capable of performing a single electron spin readout. We suggest a 
number of possible technologies that may be capable of perfomiing such a readout below. 
Before mentioning these we return to Equation (13). There we noticed that this pulse sequence 
was able to SWAP the qubit between the highest and lowest weight subspaces of the nuclear 
and electronic spins^ J = 1/2 and 5 = 3/2, and in the process increase the detectable magnetic 
signature by a factor ^ 1000. 

We will now further enquire whether there exists similar pulse sequences which could further 
SWAP out the quantum infonnation now stored in the outer electronic qubit into a separate 
coupled electronic system of larger spin. This could be very advantageous, for example^ for the 
purposes of coupling into the spin-chain, a nano-molecular magnet with a spin of 10 — 30 as a 
potential readout site [44]. As the SWAP gate is built from CNOT gates we can focus on the 
latter. Further, the pulse sequences must only consist of higjier spin representations of the Pauli 
group. There may be many such pulse sequences and we here present one set of sequences. 
The ones given below perform a CNOT between the highest and lowest weight subspaces of 
spins / and S. We set / = a/2 and S = b/2, where b>a and both a and b are odd integers. 
Letting a = [a/2] and j3 = [b/2l where [] is the integer part of the fraction, we can find 
CNOT pulses for the following four cases: 



a even, /? even: C^I, S; NOT) =^ [~~Sy][-~S, ~ I,][-^2S,I^][Sy]e-^'''^ , 
C^SJ; NOT) = [~-Iy][-^S. - I.][+25.J.][/3,]e^-/^ , 



(19) 



(18) 
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a odd, j3 odd: C^I, NOT) ^ [-Sy][~t~S, I,][+2S^I^][Sy]e'^'^/^ , (20) 
/; NOT) = hJy][+S, + I.n+^SJ^n+Iyje-^'^/^ , (21) 

a odd, ^ even: 5; ATOT) + /,][+25,/,][--5^]^^^^ , (22) 

5 q5, J; NOT) ^ [/.j^^, ^ /j[+26',J,][~^/,]e-^-/^ , (23) 

a even, ^ odd: C^I, S; NOT) :==^ [Sy][-S, - J,][-i-2SJ,][^Sy]e'^'^^^ , (24) 
J; JVOT) = + I,][+2SJ,][^Iy]e^'^^'' . (25) 

10 A Type-I quantum computer design will require the capability of single qubit readout and this 
translates in our case to the capability of reading out the spin state of the endohedral electrons. 
One can classify the various readout methods into the general categories (i) force, (ii) electric, 
and (iii) optical, measurement tecliniques. As was mentioned before, since the electronic 
wavefunction of the trapped atom is totally confined within the C60 cage, spin measurement 

15 techniques that involve the physical transport of this spin outside the molecule are not possible 
here. Some techniques, such as ODMR (optically detected magnetic resonance), have so far 
not proved possible with ^^@C60 and ^^P@C60, and thus the addition of separate readout sites 
which can be dynamically coupled and decoupled from the primary processor via NMR/ESR 
pulses are wan-anted. 

20 

Possible techniques in the above mentioned categories include (i) Magnetic Resonance Force 

Microscopy (MRFM) [47]; (iia) Micro-Squids, a technology which is already capable of 
discriminating a Arris = 30 [46]; (iib) Scanning-Tunneling-Microscope ESR, a technique 
which, though not well understood, has yielded single molecule ESR spectroscopy of iron 

25 impurities in Silicon [45], and in single BDPA molecules on Silicon [48]. (iic) A single 
molecule endohedral fullerene single electron transistor [49]. This is a very challenging 
technology which has yielded a single electron tunneling current through a C60 molecule with 
an electro-mechanical coupling to the quantized motion of the entire molecule. Performing 
such an experiment with ^^@C60 may yield endohedral spin information; (iiia) Magnetic 

30 coupling of the endohedral electronic spin to a solid-state paramagnetic optical dipole such as a 
nanocrystaline N-V center in diamond. The optical paramagnetic center can then be probed via 
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Optical shelving techniques [50]; (iiib) ODMR via a magnetic coupling to a paramagnetic 
endohedral adduct which is chemically bonded to the ^^@C60 molecule. For use with ODMR 
the adduct endohedral should ideally possesses an optical transition in the visible spectrum 
such as Bra:Sc3^JSf@C80 [51] . 

5 

The techniques suggested above are, in some cases, already at the single electron detection 
level (though not capable of detecting the spin orientation), while the out-coupling of the 
quantum information into a large-spin system may allow the techniques of micro-Squids and 
MRFM as they presently stand to act as a readout of a single endohedral electronic spin. 

10 

The direction taken in the present invention with regard to the constraction of a quantum 
information processor is to follow as closely as possible the system Hamiltonian that nature 
provides. The globally addressed approach has very significant advantages over the local 
gating approach. Global addressing avoids the buildup of numerous decoherence pathways 

15 associated with the effects of the local gates on the processor. Global addressing also avoids 
the very significant problem of the scaling up of the extemal resources required to execute local 
gating, e.g. nimierous metallic contacts or multi-fi-equency MW/RF generators with increasing 
bandwidths. Although it may be true that for local gating, the extemal resources required 
increase polynomially with the number of qubits in the processor, such an overhead may not be 

20 experimentally feasible as the size of the processor grows large. 

In some sense the technical complexities involved in the construction and operation of local 
gating are transferred into the quantum cellular automata "software". As shown above, the 
detailed operation of a quantum cellular automata quantum computer architecture does not 

25 follow the standard quantum circuit model. This should not be seen as a disadvantage as all 
quantum circuit algorithms can be "compiled" to run on a QCAQCA with a polynomial 
overhead [27]. Indeed the QCAQCA may be able to run programs that do not follow the 
quantum circuit model (see [52] for an example of a non-circuit quantum program in another 
architecture). Furthermore, as the fijndamental quantum operations are implemented via the 

30 very well developed methods of NMR and ESR, the fidelity of software execution is 

heightened. The QCAQCA also has the very important property that it can execute quantum 
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computations in parallel, a characteristic necessary for the ultimate scalability of the 
architecture. 

We have considered the endohedral fullerene material ^^@C60 and ^^P@C60 and argued that 
5 the system Hamiltonians and interactions present are sufficient to implement universal quantum 
computation via a quantum cellular automata architecture. Following the DiVincenzo criteria: 
(i) we argued that the Group-V endohedral materials behave as essentially electrostatically 
isolated nanoscopic molecular neutral atom traps that can be physically manipulated with 
relative ease using STM techniques and perhaps may be self-assembled using fullerene supra- 

10 chemistry, (ii) the electronic decoherence times are perhaps the longest seen in any molecular 
system, of the order of seconds at temperatures of ^ T'^K, which is also an upper limit for the 
electronic dephasing time. The nuclear decoherence times have not yet been measured but they 
are expected to be several orders of magnitude longer than the electronic decoherence times, 
(iii) The complete polarization of both the nuclear and electronic computational spins is 

15 feasible as half of the spins (electrons)^ within the entire system are completely polarised at 
moderate conditions of temperature and magnetic field strength. There are efficient cooling 
schemes exist to shuffle the unpolarized spins away from the computational spins, (iv) 
entanglement can be generated via the inter-molecule magnetic dipole interaction, the strength 
of which has been measured to be 50MHz. Armed with the two species ^^N@C60 and 

20 ^T@C60, we showed that one has sufficient control to implement a two-component quantum 
cellular automata quantum computing architecture storing the quantum information in the 
nuclear spins while using the electrons as a local bus. We further foimd that the expected 
"figure of merit" compares very well with other proposed soKd-state quantum computer 
designs. We also discovered pulse sequences that SWAF^d the quantum information between 

25 the highest and lowest weight spaces of two spins of different size, e.g. 5 = 1/2 ^ 5 = 11/2. 
These sequences could prove useful for out-coupling the quantum information into a spin 
readout system with large spin. 

DiVincenzo' s final criterion: (v) an efficient single qubit readout, is perhaps the most difficult 
challenge for any solid-state based quantum computer design. We have proposed various 
30 possible ensemble and single-issue (Type-I), readout technologies many of which are 

themselves the subjects of intense study. The ideas presented here combine the tremendous 
power of NMR and ESR science together with the very "clean", and almost atomic, systems 
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presented by the endohedrals ^^@C60 and ^^P@C60. These, married with supra-molecular 
chemistry 5 provide a very real physical implementation of a quantum information processor. 
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Claims: 

1 . A method of quantum computing in a system comprising a plurality of spin systems 
of two distinct types A and B arrayed in an ABABAB . . . sequence, wherein each A-type spin 
system comprises a first spin subsystem Al and a second spin subsystem A2 coupled to one 
another, and each B-type spin system comprises a first spin subsystem Bl and a second spin 
subsystem B2 coupled to one another, wherein the first spin subsystems (AI3I) of 
neighbouring A and B-t>^e spin systems are substantially isolated from one another and the 
second spin subsystons (A2,B2) of neighbouring A and B-type spin systems are coupled to one 
another, said method comprising the steps of: 

c) encoding quantum information on the first spin subsystems (A1,B1); and 

d) employing the second spin subsystems (A2,B2) as a local bus for the information encoded 
on the first spin subsystems. 

2. A method of quantum computing as claimed in claim 1, wherein the step of 
employing the second spin subsystems (A232) as a local bus comprises the steps of: 

a) transferring the information encoded in the Bl spin subsystems to the B2 spin subsystems 
via a first spin swapping operation; 

b) tagging a subset of the A2 spin subsystems^ with said tagging being effected by 
manipulating those A2 spin subsystems for which B2 spin subsystems of the nei^bouring 
B-type spin systems fijlfil a predetermined condition; 

c) transferring the infomiation encoded in the B2 spin subsystems to the Bl spin subsystems 
via a second spin swapping operation; and 

d) performing a spin manipulation on a subset of first spin subsystems (Al.Bl), with said 
manipulation being targeted to those first spin subsystems (A1,B1) for which the associated 
second spin subsystem (A2,B2) fiilfils a predetermined condition. 

3. A method as claimed in claim 2, wherein the predetermined condition of the 
associated second spin subsystem in step (d) according to which the first spin subsystems are 
targeted is the condition of being tagged or untagged following step (b). 

4. A method as claimed in claim 2 or 3, further comprising an initialisation step of 
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polarising a group of second spin subsystems to a substantially common state prior to 
performing said first spin swapping operation. 

5. A method as claimed in any of claims 2 to 4, farther comprising the step of 
returning said second spin subsystems (A25B2) to their initial configurations after performing 
spin manipulation in step (d). 

6. A method as claimed in claim 4, wherein the step of returning said second spin 
subsystems (A2,B2) to their initial configurations comprises: 

a) transferring the information encoded in the Bl spin subsystems to the B2 spin subsystems 
via a third spin swapping operation; 

b) reversing the tagging of said subset of A2 spin subsystems; and 

c) transferring the information encoded in the B2 spin subsystems back to the Bl spin 
subsystems via a fourth spin swapping operation, 

7. A method of quantum computing as claimed in any preceding claim, wherein said 
first spin subsystems (Al.Bl) are spin subsystems. 

8. A method of quantum computing as claimed in any preceding claim, wherein said 
first spin subsystems (A1,B1) are a nuclear spin subsystems IaJb- 

9. A method as claimed in any preceding claim wherein said second spin subsystems 
(A232) are spin n/2 subsystems where n is an odd integer. 

10. A method of quantum computing as claimed in any preceding claim, wherein said 
second spin subsystems (A2,B2) are electron spin subsystems Sa,Sb. 

11. A method of quantum computing as claimed in claim 1 0 when dependent on claim 
S, wherein the coupling between the nuclear and electron spin subsystems of a spin system is 
via the hyperfine interaction. 

12. A method of quantum computing as claimed in any preceding claim, when 
dependent on claim 2, wherein said spin swapping operations are hyperfine spin swapping 
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operations. 

13. A method of quantum computing as claimed in claim 10 or any claim dependent 
thereon^ wherein said electron spin subsystems Sa and Sb of neighbouring A and B-type spin 

5 systems are coupled to one another via a magnetic dipole interaction. 

14. A method as claimed in claim 2 or any claim dependent thereon, wherein said first 
spin subsystems (A1,B1) are nuclear spin subsystems IaJb having spin 1=1/2, and said second 
spin subsystems (A232) are electron spin subsystems Sa.Sb of free electrons each having 

10 S=3/2, wherein a qubit in the electronic spin subsystems is encoded as an electronic inner qubit 
defined as {|ms=±l/2}, and said nuclear spin subsystems and electron spin subsystems within a 
spin system (Ia?Sa; Ib^Sb) are coupled via the hyperfine interaction. 

15. A method as claimed in claim 14, wherein the step of employing the electron spins 
15 as a local bus for the information encoded on the nuclear spins comprises performing one or 

more h)^erfme spin swapping operations between the nuclear spins I and electron spins S by 
applying NNdR/ESR pulses according to the swap operation equation: 

20 

16. A method as claimed in claim 15, wherein a sequence of NMR/ESR pulses as 
represented in Fig. 4 is applied to the spin systems to effect said hyperfine swap operation. 

17. A method as claimed in claim 2 or any claim dependent thereon, wherein said first 
25 spin subsystems (Al.Bl) are nuclear spin subsystems IaJb having spin 1=1/2, and said second 

spin subsystems (A232) are electron spin subsystems Sa.Sb of free electrons each having 
S=3/2, wherein a qubit in the electronic spin subsystems is encoded as an electronic outer qubit 
defined as {|ms=-±3/2}, and said nuclear spin subsystems and electron spin subsystems within a 
spin system (Ia,Sa; Ib,Sb) are coupled via the hyperfine interaction. 

30 

18. A method as claimed in claim 17^ wherein the step of employing the electron spins 
as a local bus for the information encoded on the nuclear spins comprises performing one or 
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more hypei-fine spin swapping operations between the nuclear spins I and electron spins S 
according to the swap operation equation: 

19. A method as claimed in claim 1 5 or 18, wherein said swap operation is selectively 
appKed to one of said spin system types only (A or B only) by applying additional NMR/ESR 
pulses to decouple the hyperfme interaction of the other of said distinct qubit types. 

20. A method as claimed in claim 18, wherein wherein a sequnece of NMR/ESR pulses 
as represented in Fig. 5 is applied to the spin systems to decouple the hyperfine interaction of 
the other of said distinct spin system types. 

21. A method as claimed in any preceding claim, wherein each spin system is 
implemented as a fuUerene molecule having an endohedral atom or molecule. 

22. A method as claimed in claim 21, wherein said distinct types of spin system are 
obtained by providing fullerenes with two different endohedral atoms or molecules therein. 

23. A method as claimed in claim 22, wherein said different endohedral atoms or 
molecules are different atoms each having a free electron and each located in the same periodic 
group. 

24. A method as claimed in claim 23, wherein one of said spin system types is 
implemented as ^^N@C60 and the other of said spin system types is implemented as ^^P@C60. 

25. A method as claimed in any one of claims 21 to 24, wherein said fullerenes are 
positioned in said ABABAB. . . array on a substrate. 

26. A method as claimed in claim 25, wherein said substrate is silicon. 

27. A method of quantum computing in a system comprising a plurality of spin systems 
of two distinct types A and B arrayed in an ABABAB sequence, wherein each spin system 
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comprises a spin Vt. entity, said method comprising performing one or more global unitary 
operations on the A spin systems, wherein: 

Jjj is the conditional application of the unitary U when U-- NOT to the A spins 

depending on the state of A^s neighboring B spins; 
5 /is the sum of the states of the neighboring B spins^ with spin up(down) defined as 

1(0) such that [0,1,2]; 
and wherein: 

said global unitary operations are performed by applying pulse sequences according 
to one of the following sequence notations (a)-(c): 

10 

^) ^ B(B, {Xf}) . {X^i^^,} - • ^; {X^}>{^4/2} - {^;72- J ' ^; {^.^}) . 

read from right to left where 2d = d\ or 

b) M == {Z„^} - J ■ {^4/2} ■ ^; {X;:^}>{Z:^./2} • {-Y^/2^ J ' ^; {X/}) . 
15 read from right to left where 6 = 9; or 

c) implemented as an application of {^^^} to all the B spins, followed by -^2, followed by 
an application to all the B spins; 



20 wherein the controlled operation of applying t/ to an •A spin with left and right B neighbours 

such that the B spins control the unitary application is denoted by C{B, A\ U). and the function H 
is in turn defined by the equation: 

EiB,A',U)=C{Bon the left 5 the right? 

and wherein in equations (a)-(c) the NMR/ESR notation for je/2 pulses used is: 
and rotation through the angle 7 about the axis, for example, is denoted as: 

28. A method as claimed in claim 27, wherein U = 

30 

29. A method as claimed in claim 28, wherein 6 = tt. 
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30. A method as claimed in any one of claims 27-295 further comprising performing 
one or more global unitary operations on the B spin systems, said unitary operations being 
defined for the B spins depending on the state of B's neighboring A spins in terms identical to 
the unitary operations Jjj other than the substitution of the symbol B for Ji and vice versa at 
each occurrence of said symbols in each foregoing equation and definition, 

31. • A method as claimed in any preceding claim, wherein said spin V2 entities A and B 
have an Ising-like interaction defined as J j^S^ S^; ^ where Jd is a measure of the magnetic 
dipole coupling strength between the A and B entities. 

32. A quantum computer comprising: 

a) a plurality of entities providing quantum spin systems of two distinct types A and B arrayed 
in an ABABAB . . . sequence, wherein 

each A-type spin system comprises a first spin subsystem Al and a second spin 
subsystem A2 coupled to one another, and each B-type spin system comprises a first spin 
subsystem B 1 and a second spin subsystem B2 coupled to one another, 

the first spin subsystems (A1,B1) of neighbouring A and B-type spin systems are 
substantially isolated fi-om one another and the second spin subsystems (A25B2) of 
neighbouring A and B-type spin systems are coupled to one another, 

b) an NMR/ESR pulse generator associated with said array of entities for generating and 
applying pulses to said array of entities; and 

c) means for reading the spin state of one or more of said entities. 

33. A quantum computer as claimed in claim 31, wherein each spin system is 
implemented as a fiallerene molecule having an endohedral atom or molecule. 

34. A quantum computer as claimed in claim 30 or 3 1, further comprising a computer 
program associated with said NMR/ESR generator, said computer program including 
instruction which when executed are effective to cause said NMR/ESR generator to generate 
pulses effective to: 
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a) encode quantum information on the first spin subsystems (Al.Bl); and 

b) employ the second spin subsystems (A232) as a local bus for the information encoded on 
the first spin subsystems. 

5 35. A quantum computer comprising: 

a) a plurality of entities providing quantum spin systems of two distinct types A and B arrayed 
in an ABABAB. - . sequence, wherein each spin system comprises a spin i4 entity copied to 
its two neigjibours 

b) an NMRyESR pulse generator associated with said array of entities for generating and 
10 applying said pulses to said array of entities; and 

c) means for reading the spin state of one or more of said entities. 

36. A quantum computer as claimed in claim 35, further comprising a computer 

program associated with said ISDVIR/BSR generator, said computer program including 
15 instmction which when executed are effective to cause said NMR/ESR generator to gemerate 
pulses effective to perform one or more ^obal unitary operations on tihie A spin systems, 
wherein: 

A^ is the conditional application of the unitary U when U NOT to the A spins 
depending on the state of ^'s neighboring B spins; 
20 / is the sum of the states of the neighboring B spins, with spin up(down) defined as 

1(0) such that /€[0,1>23; 
and wherein: 

said global unitary operations are performed by applying pulse sequences according 
to one of the following sequence notations (a)-(c): 

25 

read from right to left where 2b = 6; or 
30 read from right to left where 6 = B; or 



SUBSTITUTE SHEET (RULE 26) 
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c) implemented as an application Spins, followed by ^2, followed by 

an application of{^'^}toallthe0 spms; 

wherein the controlled operation of applying U to an A spin with left and right B neighbours 
sixch that the B spins control the unitary application is denoted by C{B, A; U). and the function 
is in turn defined by the equation: 

E{B, ^; U) = C(Bon the left, A; U)CiBon the right 5 
and wherein in equations (a)-(c) the NMR/ESR notation for n/2 pulses used is: 

and rotation through the angle 7 about the z— axis, for example, is denoted as: 
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